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Abstract 

Making use of a simple unitary transformation we change the hamiltonian 
of a particle coupled to an one dimensional gas of bosons or fermions to a 
new form from which the many body degrees of freedom can be easily traced 
out. The effective dynamics of the particle allows us to compute its damping 
constant in terms of the reflection coefficient of the interaction potential and 
the occupation number of the environmental particles. We apply our results 
to a delta repulsive potential. 
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I. INTRODUCTION 



When one wants to study the quantum dynamics of a non-isolated system a possible way 
to attack the problem is to assume that the system we want to study is coupled to a bath 
of harmonic oscillators with a specific spectral function . When this is done we are not 
interested in the microscopic details of the interaction of the system with its environment 
but rather in those essential features that help us to recover the phenomenological equations 
of motion in the classical limit. This has been the approach used to study dissipative effects 
in the phenomena of macroscopic quantum tunneling and coherence [^J. 

A second approach to investigate the quantum dynamics of a system interacting with its 
environment would be to treat the problem from a microscopic point of view. In this case 
there are no standard expressions for the effective dynamics of the system of interest and 
for each problem they must be developed from the microscopic equations according to the 
appropriate approximations for the particular interaction we are dealing with. There are 
very important problems which are suitable for the microscopic approach, such as; motion of 
ions in normal [|]|| and superfluid He 3 [|| , mobility of particles in metals and diffusion 
of heavy particles in solids [[|, to name just a few. 

Actually we could also incorporate in the microscopic approach the recently developed 
formulation for dealing with mobility and diffusion of quantum solitons in problems of cou- 
pled |]10| or non-linear field theories [|TTJ. Here, the collective coordinate method fll"2 



is 



shown to lead us to another very simple model |13[ for describing quantum dissipation. 

In this paper we have shown that the hamiltonian of a particle coupled to bosons and 
fermions in one dimension can also be mapped into the model of ref. ||13|| . In this case it is 
simply done through the application of a unitary operator which uniformly translates the 
positions of all the environmental particles changing the origin of their coordinate frame 
to the position of the external particle. As a matter of fact this procedure is valid for any 
interaction potential which depends only on the distance between the particle of interest 
and each of the environmental particles. At the end we specialize on the delta repulsive 
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potential because of mathematical convenience and transparence to analyze the effects of 
foward and backward scattering processes. 

We shall start by applying the above-mentioned unitary transformation to the origi- 
nal hamiltonian and showing how its transformed version can be used in the well-known 
Feynman- Vernon approach |I],|nj for dealing with the reduced density operator of the exter- 
nal particle. This is done in section II. 

In section III we briefly review the coherent state representation for the density operator 
of the system and obtain the real time reduced density operator of the particle from which 



we can get its damping constant. Although this has been presented in another paper [10 
we decided to include it here to make this contribution more self-contained and to correct 
some of our previous expressions. 

In section IV we analyze the expression for the damping constant and obtain its temper- 
ature dependence in the case of bosonic and fermionic environments and finally present our 
conclusions in section V. We have also included an appendix where we discuss the approach 
in thermal equilibrium and the approximations made in our calculations. 

II. THE CANONICAL TRANSFORMATION 



Let us start by describing the model that will be used throughout this paper. We assume 
that a particle of mass M is interacting with an one dimensional gas of N fermions or bosons 
of mass m through the hamiltonian 

2=1 2=1 

where p, q, pi, and are respectively the momenta and positions of the heavy particle and 
of the i th environmental particle and U(q — qi) can, in principle, be any interaction potential. 
We shall only assume it has a delta repulsive form later on in this work. 

The dynamics of the particle is governed by its reduced density operator which is given 

by 0,0 



3 



p(x,y,t) = J dq 1 dq 2 ...dq N (x,q 1 , ...,q N |e*»*p(°) e ^ \VAu -An) (2.2) 

where p(0) is the initial density operator of the whole system. 
Now, introducing the unitary operator U 

U = e~*L'i=i m (2.3) 

we can use the identity U~ X U = 1 between any two operators in (2.2) to get 

p~(x,y,t) = [ dq 1 dq 2 ...dq N (x,q 1 ,...,q N \U~ 1 e lI 7rp'(0)e^ I TrU\y,qi,...,q N ) (2.4) 



where H' is the transformed hamiltonian [15 



and p'(0) = U p(0)U~ l is the transformed initial condition which we analyze below. 

It should be noticed that the application of the unitary transformation U to the hamilto- 
nian H has modified it in such a way that now the external particle interacts via momentum- 
momentum coupling with a set of particles which are scattered by a potential U(q). The 



same transformation was used by Castella and Zotos |L6| for the exact diagonalization of 
the problem via the Bethe ansatz. 
Now, inserting the closure relation 

J dxdq 1 dq 2 ...dq N \x,q 1 ,...,q N )(x,q 1 ,...,q N \ = 1 

twice between the unitary operators U and U~ l and the transformed density operator p'(t) 
one has 

p{x,y,t) = J dq 1 dq 2 ...dq N dr 1 dr 2 ...dr N (r x , ...,r A r|e"' r ^= lPj{:!: " s/) |gi, ...,q N ) 

(x,qi, ...,q N \p'(t)\y,r u ...,r N ) (2.6) 

which is more cumbersome than (2.2) due to the existence of the first term on its r.h.s.. 
This complication was brought about by the unitary transformation and reflects the fact 



that after solving the dynamics in the more appropriate representation we need to return to 
the original representation from where we started. 

Nevertheless things are not so awkward as they seem. Since we are most of the time, 
and particularly in this work, interested in the average values of observables such as (x(t)), 
(p(t)), (x 2 (t)) or (p 2 (t)) we do not really need (2.6) as it stands but rather its diagonal 
elements only. 

The justification for this assertion for observables which are only x-dependent is trivial 
because 

(x n (t)} = J dxx n p(x,x,t) (2.7) 

and the first term on the r.h.s. of (2.6) becomes 5(qi — ri)...5(qN — r N ) which can be trivially 
integrated yielding 

(x n (t)} = J dxdq 1 ...dq N x n (x,q 1 ,...,q N \p'(t)\x,q 1 ,...,q N ). (2.8) 

For p-dependent operators the demonstration is also easy although not immediate as in 
(2.7-2.8). Now, what we need is 

«'» = i * n — s (2 ' 9) 

which upon the application of the derivatives to (2.6) yields 
(p(t)) = ^Lim x ^ y J dq l dq 2 ..Aq N dridr 2 ..Ar N {r u ...,r N \e~^^ Pj{x ~ y) \q 1} ...,q N ) 

In the limit x — > y the first term of the integrand above becomes once again 8{q\— ri)...S(qN — 
r N ) which can be integrated to give 

(p(t)) = ^:Lim x ^ y J dq x dq 2 ...dq N ^ ~ ^ ( x > 9i> •••> QN\p\t)\y, qi, q N ) (2.11) 

and so, what we really need in both cases is, 

P'(x,y,t)= J dq 1 dq 2 ...dq N (x,q 1 ,...,q N \p'(t)\y,q 1 , ...,q N ). (2.12) 
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Here one should notice that this result is valid only if we want to compute the above- 
mentioned averages. Were we interested in the time evolution of p(x, y, t) this statement 
would not be valid. 

Before leaving this section, let us try to clarify the choice we will make for the initial 
condition p'(0) in this work. As we have seen from (2.4) p'(0) can be prepared in terms 
of the eigenstates of the transformed hamiltonian (2.5). Our particular choice will be to 
assume that the particle is in a pure state (a wave packet centered at the origin) and the 
environment is in equilibrium in the presence of the particle. This means we can write p'(0) 
as 

P'(0) = PsP' eq = p s e- pHe (2.13) 
where p s refers only to the particle and H e is given by 

tf« = g(|| + t%)) (2-14) 

Notice that this is only valid as long as one neglects the spread of the initial wave packet 
about the origin when preparing the initial state of the composite system. 



III. THE FEYNMAN- VERNON APPROACH 

Having specified the central object of this paper we should now choose the most ap- 
propriate representation in which we should write it. Since our environment is composed 
of indistinguishable particles (bosons or fermions) it is convenient to write the transformed 



hamiltonian (2.5) in its second quantized version which reads [17| 



where = h (i\p'\j) is the matrix element of the momentum operator of a single environment 
particle between eingenstates \i) and \j) of if e in (2.14) , fij are the eigenfrequencies of these 
states and p is the chemical potential. The operators and a\ are the standard annihilation 
and creation operators for bosons or fermions. 



Then we can use the initial conditions (2.13) in (2.12) to write the standard Feynman- 
Vernon expression fliyi4]| 

p(x, y,t) = J dx J dy'J(x, y, t; x', y', 0)p s (x', y', 0) (3.2) 

where 

J(x, y, t; x', y', 0) = J dq 1 ...dq N dr 1 ...dr N ds 1 ..ds N (x, q 1} q N \e~% H,t \x' , r u ...,r N ) 

(r x , r N \e~ pHe \s x , s N )(y', s 1 , s N \e^ H>t \y, q u q N ). (3.3) 

However, as we have chosen to represent our hamiltonian H' in its second quantized form 
we had better perform all the integrals in (3.3) in the coherent state representation for the 
creation and annihilation operators |0| involved therein. In doing so we can rewrite J in 
(3.3) as 

J(x, y ,t;x',y',0) = I — ^ — (x,a|e-^'V,7>(7|e- /JH »|^(y',5|e^|y,a>. (3.4) 

J IX TX TX 

where a is a complex vector with infinte number of components and — = 
- Hi d(Reoti)d(Ima[i). Equation (3.4) can be written in a more appropriate form if we 
use the functional integral representation for all the time evolution operators (in real and 
imaginary time) of its integrand. We can write 



J(x, y, t; x', y', 0) = / Vx(t')Vy(t')e^ s °^~ Solv] ^[x, y] (3.5) 

where 

S [x] = f dt'\Mx 2 {t') (3.6) 
Jo 2 

is the action of the particle if it were not coupled to the environment, J-[x, y], the well-known 
influence functional, is given by 

d 2 ad 2 ^d 2 S ? /•«* 



F[x,y] = f—-^— PR (rJ) I V 2 a />? e K*[^-^) (3 . 7 ) 

J IX IX 7X Jy J 8* 



with 



Si[x,a)= f dt' \-^^2(a* n a n - a n a* n ) + xJ2 n 9m, n a* m a n - ^2(hn n -/x)a*a„] (3.8) 

JO \ £ „ r „ n n I 



and 

exp {e-M Kl »-ti'ft6n} 



Z 

d 2 a 



Z = j ^exp{e-^ nn ^\a n \ 2 }. (3.9) 

The integral in (3.7) is exactly the same that has been evaluated in the context of 
the polaron dynamics in ref. [|L0[| . The only difference is that we now must deal with the 
environmental particles as being massive bosons or fermions. The result of the integration 
follows exactly as in |H| and reads [|19[ 



r[x,y] = [det(lTNr[x,y])F 1 (3.10) 

where the —1 refers to bosons and +1 refers to fermions, N is the ocupation number operator 
given by 

Nij = <),j7T; 

* = eft*-*) zp 1 ( 3 - U ) 
and the functional F[x, y] is such that 

T nm [x, y] = W nm [x] + W* nm [y] + £ W* nk [y]W km [A (3-12) 

k 

with W nm satisfying the integral equation 

W nm [x, t] = f T dt' W° m (x, t') +J2 f T dt ' f dt " W nk(i, t')W km (x, t") (3.13) 
Jo k Jo Jo 

where 

Km(t') = ignmXit'y^-^ '(1 - S nm ). (3.14) 

So far our results are exact. 
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Since the influence functional T must still be integrated in (3.5) its form (3.10) is useless. 
In order to evaluate the double functional integral in (3.5) we have to make some approxi- 
mations in our present expression for T\x, y\. Observing, from (3.13) and (3.14), that T[x, y] 
depends on the velocities x and y we can rely on the assumption that only slow paths will 
contribute to (3.5) and expand (3.10) up to quadratic terms in x and y. 

This was the procedure followed in reference 0] and therefore we shall only quote the 
final result of this approximation which reads 

j=\x, y\ = e -i^,y] e -k^n[x,y] ( 3 _ 15 ) 

where 

f dt"(x(t') -yit'^hT^t' -t n )(x{t v ) + y{t n )) (3.16) 
jo Jo 

and 

®r= fdt' f dt»{x{t')-m)W R (t'-t"){x{tn-y{t v )) (3.17) 

JO JO 

with 

r i?(<0 = 2 Yl \9ij\ 2 (ni + n j + ZniUj) cos(Qj - Qj)t (3.18) 

and 



r /00 = H \9ij\ 2 ( n i ~ n j) sin(fii - Qj)t. (3.19) 
From these expressions we can determine the damping function j(t) as 

7(t) = ~M~tt = ~2M J du J ^' S ^^')H^)-n(^')\{^-^')cos{u-u } ')t (3.20) 

and the diffusion function D(t) as 
d?Y h 2 f f 

D(t) = -h-^T = YJ duJ J duJ ' S ( u > u ') H w ) + + M^M^O] (oo - Jf cos(w - u')t. 

(3.21) 
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In the expressions for ^(t) and D(t) we have introduced the scattering function S(u,u') 
which we define as 

S(w, oj') = \9a\ 2 ^ - ni)5(w' - Qj). (3.22) 
In this paper we shall be mostly concerned with the damping function (3.20). 

IV. THE DAMPING CONSTANT 

Having deduced the explicit form of ^(t) our next step is to investigate its long time 
behavior. As in ref. JTOf it is our goal now to evaluate the scattering function S(u,u') for 
a specific choice of our model hamiltonian. Let us assume that the interaction potential 
between the external body and the environmental particles is given by 

U(q - ft) = V a6(q - ft) (4.1) 

where Vo and a are, respectively, the height and the range of the potential barrier. 

Now, because of the unitary transformation we used, the expression of S(u,u f ) depends 
only on the evaluation of the matrix elements of the single particle momentum operator 
between two different states in (2.14). As this is a hamiltonian of non-interacting particles 
subject to a delta repulsive potential at the origin it is a very simple task to obtain its eigen- 
states. Observe that this problem is essentially the one of an impurity in one dimension that 
has been studied using many different techniques, from bosonization and renormalization 



group calculations P0| , |2l|| to Bethe ansatz approaches [22]. 

Let us suppose that our environmental particles are contained in a very large one dimen- 
sional box of length 2L. It is then easy to show that we have even and odd solutions given 
by 0, 

UE(k, x) = — ^= cos(A;|x| + 5 k ) (4.2) 

uo(k, x) = — ■=. sin(foc) (4.3) 
yL 
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where 



with 



1 , 

-tan 
z 


( im(t + 


1 fiW/ 4- 


t\ k ) = 


IK 


2\k\D(k) 


*(*) = 


l + r(k) 


£>(*) = 


IK 

2\k\ 


k = 


4mV 


h 2 ' 



(4.4) 



(4.5) 

Notice that the phase shifts of the odd solutions are zero and that we can further relate 5 k 
to the transmission and reflection coefficients of the delta repulsive potential as (the case of 
a general potential without bound states is treated in the Appendix) 



T(k) = \t{k)\ 2 = cos 2 4 
R(k) = \r{k)\ 2 = sin 2 4- 

Our first step to get an analytic expression for S(u,u') is to evaluate 

9kk> = J L -j^ u E{k, x) ^-i^-uo(k', x)"j 

which gives 



(4.6) 



(4.7) 



i f kk' 



9kk' — — j 
+ 



L [k + k' 
kk' 



sin(k - k')L r (1 - cos(A; - k')L) 

cos 5 k ; - sin Ok 



k-k' 



k-k> 



k-k' 



sin(k + k')L . (1 - cos(A; + k')L) . . 

COS Ok ; — sin Ok 



k + k' k + k' 

When we take the limit L — > oo in (4.8) and use that, since g is purely imaginary, g^k 
~9k'k, we have 

sin 5k + sin Sk 



(4.8) 



gkk' 



— I — 5{k 2 — k 2 )kk'[cos 5k — cos 5k'} — V 



(4.9) 



k 2 - k' 2 

with V meaning the principal value. This expression must now be used in the continuum 
version of (3.22) that reads 
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7-2 

S(ui, u') = — I dk I dk'\g kk ,\ 2 5(u; - Q k )5{uj' - Q k ,) (4.10) 



7T 2 



where = Evaluating (4.10) one gets 



2m 

(4.11) 



sin 5(u) — sm6(u') 



uj — uj' 



Now we can finally evaluate the integral for j(t) in (3.20) which in the long time approx- 
imation gives us 

7 (0 = 7CW) ( 4 - 12 ) 

where 

is the damping constant that we want to study. In this expression the reflection coefficient 
R(E) can be easily obtained from (4.5) as, 

m = ^ (4.i4) 

where E = " is what we shall be calling the strength of the delta repulsive potential. 
Let us now study (4.13) for fermions and bosons at zero and finite temperatures. 

A. Fermions 

At zero temperature it is straightfoward to write 7 for the fermionic case. As we know, 
at T = 0, 

£L = -6(E-E F ) (4.15) 
where E F stands for the Fermi energy of the system. Therefore, 

-(T = Q) = — E o E f (Am) 

Equation (4.16) tell us that for strong interactions or low densities (Ep « E ) the 
damping constant is dominated by E F . Although the reflection coefficent tends to one at 
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very low energies, the number of fermions being scattered tends to zero. On the other hand, 
for weak interactions or high densities (Ep » Eq) the expression for the damping constant 
is governed by Eq because at very high energies the reflection coefficient is vanishingly small 
and, consequently, all the scattering effects come from the low energy fermions (E < Eq). 

The finite temperature behavior of (4.13) will now be studied in two limits; the low and 
high temperature regimes. For studying the low temperature limit we make use of the well 
known Sommerfeld expansion |24| which allows us to write 



1J ' vrv 12 dfi 
where, for the one dimensional Fermi gas, 

n 2 (k B T) 2 

Using (4.14) and (4.17) in (4.13) and keeping only terms 0(T 2 ) we have 

~(T) - mE ° X Ji _i_ 71-2 1 ~ x ( k zT x ~ 
1[ '~ 27rhMl + x\ 12(l + x) 2 \E F 

x = -f (4.19) 
Eq 

From (4.19) we see that the finite temperature correction to 7 is strongly dependent on 
x. For strong interactions or low densities (x < 1) the low temperature correction to 7 is 
positive. This means that we are starting to populate states with E > Ep in a region where 
the reflection coefficient is still appreciable. 

For weak potentials or high densities (x > 1) we have already seen that the damping 
is dominated by low energy states. As we increase the temperature there is no change in 
the behavior of the high energy fermions because of the very small values of the reflection 
coefficient in this region and only the depletion of the low energy states will influence the 
behavior of the damping constant lowering its zero temperature value. 

The classical limit, which we define as k B T » Ep, can also be analyzed if we use that 
in this limit the Fermi occupation number is the same as in the distinguishable particle case, 
namely ||17|| ; 
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n(E) = ±J™?l e -f>E (4 . 20) 
gym 

where g is the degeneracy factor and A the linear density of fermions, 

(, 21 , 

Using (4.20) and (4.13) we can approximately evaluate that integral in two limits of 
interest, 



. m k B TE F ^ , _ 



^ = 1FS5» E ° <<k ° T - (422) 

Actually there are three distinct cases which we can analyze in the whole temperature range. 

a) E F » E ; in this case 7(0) is proportional to E and starts to decrease as the 
temperature is raised. It keeps on decreasing up to high temperatures (k B T » E F >> E ) 
where it has a T -1 / 2 behavior. 

b) Eq » Ep] now 7(0) is proportional to Ep and increases with temperature. However, 
there must be a maximum value for 7(T) since its predicted behavior at high temperatures 
(k B T » E » E F ) also obeys the T" 1 / 2 law. 

c) Impenetrable external potential, E — > 00; at low temperatures the behavior of the 
damping constant is exactly the same as the previous case. However, as the temperature 
increases it will never get to the T _1//2 behavior because E >> k B T >> E F . Therefore it 
presents a T 1 / 2 for high temperatures. 

The T 1 / 2 behavior of a) and b) above is easily understood because in both cases the 
energy k B T exceeds the strength E of the potential barrier causing a strong depletion of 
the states which are sensitive to the presence of the scattering center. 

For the impenetrable potential the reflection coefficient is always one and therefore 7(T) 
is a monotonically increasing function of T and behaves as T 1//2 in the classical limit. These 
temperature dependences are illustrated in Fig.(l) and Fig. (2). 
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Ep E q % ^ 

FIG. 2. Damping constant j(T) as a function of temperature for Eq > Ep 
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B. Bosons 



The same analysis can be carried on in the case of massive bosons. The steps are not so 
straighforward as in the fermionic case but they present no special difficulty either. It is a 
very easy task to show that the occupation number for bosons has the following assymptotic 
limits; 

n(E) = - B , , E + « k B T 
ti + \n\ 

n(E) = E + » k B T. (4.23) 



Observe that for bosons the chemical potential /i is always negative JT? 
Now, if we assume that, when T — > 0, /i(T) behaves as p5f 



A*CO = (4.24) 

we can deduce from (4.23) that 

n(E = 0) = N (4.25) 

where N ist the total number of bosons in the system. 

In order to further investigate the crossover from the low to high energy behavior of n(E) 
in (4.23-4.24) let us assume that we are in the high temperature regime which means that 



we can replace the chemical potential by its one dimensional classical expression |L7 



^(T)^^ln(|) (4.26) 
where Tq = 2r ^^ 2 is the temperature to localize a particle in a region of length 2L/N and 



A is, as in (4.21), the linear density of particles. Notice that by high temperatures we mean 
T > T Q . In this case the crossover takes place at the energy 



E c = k B T - \fj,\ = k B T 



i + V- 

2 VT 



(4.27) 



Since ln(To/T) is a monotonically decreasing function of T we see that as we increase T the 
crossover energy E c is pulled down to lower values and eventually reaches zero. This means 
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that at this temperature the occupation number can be treated classically for any energy 
E. It is a trivial matter to calculate this temperature because when is is reached one has 
E C (T C ) = or, 

T c w e 2 T . (4.28) 

Therefore, this is the temperature above which the system is classical. 

The same analysis tells us that when T « T the crossover energy E c is very high 
and, consequently, the low energy form of n(E) dominates any physical phenomenon in this 
temperature range. Actually we can use (4.23) to compute N (T), the number of particles 
in the lowest energy state, for very low T. We can do this by writing that the total number 
of bosons can be approximated by the expression 



L 2m fk B T-M dE k B T L 2m k B T 
iV ~ — \ —rr / — ^ — r-r = — \ — ^tan 



7T V h 2 Jo 2y/EE+\ij\ it V h 2 \ 

and then, assuming that 



k B 7 






fi 









(4.29) 



LimN^^T-tO-TjT = °> ( 4 - 30 ) 
which is in agreement with (4.24), one gets 

MT) « (4.31, 

where JV„(T) « ^. 

This is the result that enable us to compute the low temperature behavior of 7(T). This 
is done by integrating (4.13) by parts and using (4.23) to write 

m = ^_ /'" T -"' 1 dE ( _£2_) 2 JSL. ,4.32) 

' v ; 2nhM Jo \E + eJ E+\u\ K J 



which with the help of (4.31) yields 



From this formula we see that the damping constant vanishes at T = but presents a 
non- analytic behavior as a function of temperature in this limit. This is due to the fact 
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that n(E) is proportional to VES(E) when T — > and therefore the weight of the states 
influenced by the external particle vanishes in this limit. As the temperature is raised there 
is a fast increase of this weight. 

The high temperature limit of the damping constant has the same form as for the 
fermionic case (as it should) with the only difference that E F must now be replaced by 
ksTo. Then, 

Once again following our procedure in the fermionic case we can also study distinct temper- 
ature dependences of the damping constant as we show below. 

a) E » ksT; here 7 vanishes at T = and increases with temperature as in (4.33). 
When ksT » E >> ksT it decays as T" 1 / 2 and, consequently^ must present a maximum 
at intermediate temperatures (E » ksT » A;_bT ). 

b) E — > 00; the low temperature behaviour is the same as before but now one can never 
go beyond the regime E >> ksT >> ksT and therefore 7 behaves like T 1 / 2 . 

Notice that we have not analyzed above the condition /cbT >> E because it does not 
make sense for bosons. In this case the strength of the barrier is below the minumum energy 
required to localize a boson within 2L/N and so this would correspond to an undamped 
motion. 

The interpretation of the high temperature behavior of 7 follows similarly to what we 
have said for fermions. 

Here we can also scketch the behavior of 7 for a fixed T as we show in Fig. (3). 
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yii) 




K B T E„ 

FIG. 3. Damping constant j(T) as a function of temperature for UbTq < Eq 

V. CONCLUSIONS 

In this paper we succeeded in writing down an expression for the damping constant of a 
particle coupled via a delta repulsive potential to non-interacting bosons and fermions. 

We were able to show that using a very simple unitary transformation we could map the 
original hamiltonian into a new one where the external particle interacts via a momentum- 
momentum coupling to a new system where the particles are acted by an external potential. 
Moreover, we could also show that within the Feynman- Vernon picture all we needed was 
to find the time evolution of the reduced density operator for the transformed system. 

Although our results were quite general very little can be done beyond the gaussian 
approximation of the influence functional. On the other hand this is exactly the form we 
should have got had we started with a bath of harmonic oscillators [l],[J and therefore, a 
question still remains; why haven't we started with this simpler model? 

The answer is simple. In our present problem our input is the microscopic parameters 
of the interaction while in |l],|2| all one has is a constraint to recover the Langevin equation 
in the classical limit. There, the microscopic constants all combine in order to establish the 
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low frequency behavior of the spectral function J(uj) |T],|2| to be the form rju, where 77 is the 
damping constant, whereas in the present approach the damping constant is obtained as a 
function of the microscopic parameters and the temperature. It is clear one could try to find 
a temperature dependent spectral function J(u,T) in order to recover our function j(T). 
However, it is hard to imagine that this could be done a priori, without the knowledge 
of the specific form of j(T). Even if we can phenomenologically describe j(T) in a vast 
temperature range it is many times of great importance to understand this behavior from 
first principles. 

Notice that in the present formulation we have also a prescription to go beyond the 
gaussian approximation. This can be used if one wants to test the stability of the "mean 
field" solution to this problem. This project is being carried forward by one of us in a 
separate contribution f2~6fl . 
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APPENDIX A: FORMULATION IN THERMAL EQUILIBRIUM 



In this appendix we extend the real time calculation of the paper to imaginary time and 
stress the approximations used in the bulk of the paper. 

We look here at the partition function for the initial problem which can be written as, 

Z = tr (e~? H ) = tr (e^ H ') (Al) 

where we have used the unitary transformation (2.3) and the fact that the trace is invariant 
under this transformation. 

The partition function can be written in terms of path integrals for the external particle 
and the fields of the bath as 

Z = J DqJ D^*D^ exp (S E [q, *]) (A2) 

where the action associated with the hamiltonian (2.5) reads (h — 1), 

Se= fdtr-^-- ix\2g kk/ %V k/ - V I ^ - E k 1 %* k ) ■ (A3) 



o 



Since we are interested only in the dynamics of the particle we can trace the environ- 
mental particles exactly in order to get, 

Z = J DqeM-Seff[q}} (A4) 

where the effective action for the particle is non-linear and reads 

ft 3 Mx 2 

Seff= / dt— =F In det W[q] (A5) 
Jo I 

where the minus sign is for fermions and the plus sign for bosons. The operator W is written 
in terms of its matrix elements as, 

d 



(kt\W\k't') = 5(t-t') 



^ E k j 5 k y - ixg kk i 



(A6) 



Clearly the problem (A5) is extremely complicated and is equivalent to the one obtained 
in real time from expression (3.10). However we are interested only in the case where the 
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velocity of the particle is small or its mass is much larger than the mass of the environmental 
particles. Therefore, we look at the saddle point equations for the action (A5) in the limit 
where M — > oo. It is very easy to show that within the mean field approximation the saddle 
point equations read 

i ^ ( d 



*± s g 9 w(^G(M;^)j^ = o (A7) 



and 

( d 



h E k j G(k, t; k', t') =f ix(t) £ 9kk»G(k", t; k', t') = 6 kk ,6(t - t') (A8) 



In the limit M — > oo we choose the particular solution 

x(t) = 



j t - E k \ G (k, t; k', t') = 5 kk ,5(t - t') (A9) 



which can be easily solved by 

G (k,t;k',t') = d -^Y,h V (A10) 

P n lUJ n — &k 

where u n = for bosons and cu n = ' K C 2n + r > for fermions. 

Now we look at gaussian fluctuations around this mean field by expanding the determi- 
nant in (A5) as 

lndet W = lndet(# + #/) ~ hidetG^ 1 + tr(G #j) - -tr(G #/) 2 (All) 

where HqGq = 1 and Hj is the second term on the r.h.s. of (A6). 
Thus, the effective action reads 

rP Mr 2 rP rP 
S eff ~ S°ef f + / + dt dt'F(t - t')x(t)x(t') (A12) 

J JO JO 

where 

F(t) = J2\9kk>\ 2 G (k,t)G (k',-t). (A13) 

k,k' 
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It is useful to consider the action (A12) in Matsubara representation [[17] 

A5 e// = 4 E ( M ^n + ^K)) * W^-Wn) (A14) 
P n 

where 

k,k' lL ° n + -^fc — ^fc' 

The actual expression for F{ui n ) now depends on the form of the coupling constant gkk'- 
In the paper we considered the case of a delta repulsive interaction but in the following 
subsection we will consider the case of a general potential without bound states. This 
approach generalizes the results of section IV. 

1. Matrix element 

Let us consider the case where the potential of interaction between the particle and 
the environmental particles is localized, that is, V(x) = for x > \a\. In this case the 
Schrodinger equation for the impurity problem in the assymptotic limit can be written as 

^ k ( x ) + k 2 ^ k (x) = (A16) 



y k 2m ' 



where 

The general solution of this problem can be written in terms of even and odd wavefunc- 



tions as follows [23 



*fc ( x ) = y \ cos ( k \ x \ + <M*0) 

*£(ar) = ]J^rsgn(x) sin (k\x\ + 5 (k)) (A17) 

where sgn(x) = +1(— 1) if x > 0(x < 0), and Sg^k) and So{k) are the even and odd phase 
shifts, respectively. 

The phase shifts are related to the reflection and transmission coefficients by a simple 
relation 

23 



e 2i5 E (k) _ e 2iS (k) 

r ( k ) = g 

p 2i5 E (k) , 2i6 (k) 

t(k) = | , (A18) 

which can be inverted in order to give 



= 2 tan 



.! flm(t(k)+r(k)) 
\Re(t(k)+r(k)) 



The reflection and transmission coefficients can be now expressed in terms of the phase shifts 

as 

R(k) = \r{k)\ 2 = sin 2 (S E (k) - 5 {k)) 

T{k) = \t{k)\ 2 = cos 2 (8 E (k) - 5 {k)) . (A20) 

Thus, any information from the scattering process is encapsulated in the phase shifts. 

From (3.1) it is easy to sec that the coupling constant gkw only couples states with 
different parities, that is, 

9kk' = 9kk' = 

0,E E,0 r\ni\ 

9kv = -9k'k ( A21 ) 

The calculation of the matrix element follow the same steps as in (4.7) and (4.8) and reads 

o,e_ * kk' ^ sin (5 (k) - 5 E (k')) + sin (5 (k>) - 6 E (k)) \ 
9kk ' ~ 2Lk^k* V \ k~k> ) ■ [A22} 

The integral we need to evaluate in (A15) can now be written as 

, . _ 1 y> (sin(fl (fc) -5 E (k')) + sm(5 (k') - 5 E (k))) 2 (n k - n k >) 

[UJn) (f - (k'mEk - E k , + kj n ) ■ [ 6) 

This integral can be done by choosing a convenient countor in the complex plane which 
avoids the double pole. The final result is 

^») = r^i:^w(-f)- (a*) 
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Thus, from (A24) and (A14), the effective action reads 

= ~a ( u l + T^nl) x(u n )x(-u n ) (A25) 

P n 

which is the usual action for the motion of a particle in a dissipative environment |2| with 
a damping constant 

7=^& 2 *M(-f). (A26) 
Equation (A26) can be rewritten using the usual definition of the density of states 

D{E) = W 8 (E-E P ) (A27) 
L p 

as 

^^L dED ^ ER ^(- d T-) (A28) 

where p 2 = 2mE and v (E) = It turns out that in one dimension the density of states is 
simply 

D{E) = 2^E) (A29) 
and (A30) becomes identical to (4.13). 
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